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ABSTRACT 


4)' 

The  study  presents  a  new  integral  identity  for  the  velocity 
potential  of  three-dimensional  flow  about  a  ship  moving  with 
constant  speed  in  regular  waves.  This  integral  identity  is 
valid  outside,  inside,  and  exactly  on  the  surface  of  the  ship, 
and  is  equivalent  to  the  set  of  three  classical  identities  valid 
strictly  outside,  inside,  and  on  the  ship's  surface,  respectively. 

For  the  usual  problem  of  ship  motions  in  a  regular  sea,  the 
integral  identity  obtained  in  this  study  yields  an  integro- 
differential  equation  for  determining  the  velocity  potential  on 
the  ship's  surface,  A  recurrence  relation  for  solving  the  pro¬ 
posed  new  integro-dif f erential  equation  is  presented.^ 

N 
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1.  INTRODUCTION 

Motions  of  a  ship  advancing  with  constant  average  speed  in  regular  waves  are 
predicted  theoretically  by  using  approximate  theories  based  on  the  slenderness  of 
ship  forms.  These  theories  are  the  strip  theory,  most  useful  in  the  short-wavelength 
regime,  and  the  complementary  low-frequency  slender-body  theory;  these  complementary 
slender-ship  approximations  have  recently  been  united  and  extended  in  a  unified 
slender-ship  theory  valid  for  all  frequencies.  A  detailed  mathematical  presentation 
and  historical  account,  including  extensive  references  to  the  relevant  literature, 
of  these  slender-ship  theories  may  be  found  in  Newman.^" 

Agreement  between  strip-theory  predictions  and  experimental  measurements  has 
been  found,  in  a  large  number  of  cases,  to  be  sufficient  for  many  practical  purposes. 
This,  and  the  relative  mathematical  and  computational  simplicity  of  the  strip  theory, 
have  made  that  theory  the  most  widely  used  method  for  predicting  ship  motions. 

Indeed,  the  theory,  with  the  improvements  of  the  recently  proposed  unified  slender- 
ship  theory,  seems  likely  to  continue  to  provide  a  very  useful  and  practical  tool 
in  the  future,  even  if  significant  improvements  in  computer  performance  are  made  and 
calculations  based  on  a  three-dimensional  theory  become  more  practical. 

Notwithstanding  its  many  merits,  the  strip  theory  evidently  has  limitation,  and 
in  some  cases  there  is  a  need  for  potentially  more  accurate  calculations  based  on  a 
fully  three-dimensional  theory.  For  instance,  three-dimensional  calculations  would 
be  useful  for  predicting  the  pressure  distribution  on  a  complex  bow  shape,  such  as 
one  equipped  with  a  bulb  or  a  sonar  dome. 

At  present,  fully  three-dimensional  calculations  represent  a  difficult  task, 

and  indeed  only  a  very  limited  number  of  numerical  results  have  been  obtained  by 
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only  a  few  authors:  Chang,  Guevel  and  Bougis,  and  Inglis  and  Price.  Further¬ 
more,  these  sets  of  numerical  results  are  not  entirely  consistent:  while  agree¬ 
ment  is  good  for  some  hydrodynamic  coefficients,  discrepancies  are  very  large  for 
other  coefficients.  This  lack  of  consistency  suggests  that  the  accuracy  of  three- 
dimensional  calculations  may  be  difficult  to  control,  as  is  the  case  for  the  problem 
of  wave  resistance. 

Three-dimensional  calculation  methods  are  based  either  on  numerical  solution  of 
an  integral  equation  for  the  velocity  potential  or  on  related  assumed  distributions 
of  singularities  (sources  or/and  dipoles)  on  the  ship  surface.  Different  integral 

*A  complete  listing  of  references  is  given  on  page  27. 
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equations  can  thus  be  formulated,  and  these  can  be  solved  in  several  ways,  in  parti¬ 
cular  by  using  an  iterative  solution  procedure  or  inverting  a  matrix  of  influence 
coefficients.  The  performance  of  a  three-dimensional  calculation  method  (measured  in 
terms  of  accuracy  control,  computing  times,  and  complexity  of  implementation)  must 
obviously  depend,  to  a  large  extent,  upon  the  form  and  mathematical  properties  of 
the  integral  equation  and  upon  the  solution  procedure  selected  as  the  basis  of  the 
calculation  method.  It  thus  may  be  useful  to  consider  various  alternative  integral 
equations  and  solution  procedures. 

The  object  of  this  study  is  to  present  a  new  integro-dif f erential  equation  and 

a  related  recurrence  relation  for  determining  the  velocity  potential.  The  results 

given  in  this  study  generalize  those  obtained  previously  for  the  particular  problems 
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of  wave  radiation  and  diffraction  at  zero  forward  speed,  ship  wave  resistance, 
and  potential  flow  about  a  body  in  an  unbounded  fluid. ^ 

The  integro-dif f erential  equation,  defined  by  Equations  (5.1)-(5.4),  is  an  equa¬ 
tion  for  the  velocity  potential  <J>,  rather  than  for  the  density  of  a  related  distri¬ 
bution  of  sources  or  dipoles.  This  equation  involves  both  a  waterline  integral 
(i.e.,  a  line  integral  around  the  intersection  curve  between  the  mean  hull  surface 
and  the  mean  sea  plane),  as  in  the  problem  of  ship  wave  resistance, 10  and  a  water- 
plane  integral  (i.e.,  an  integral  of  the  Green  function  over  the  portion  of  the  mean 
sea  plane  inside  the  mean  hull  surface),  as  in  the  problem  of  radiation  and  diffrac- 
tion  at  zero  forward  speed.  The  highly  singular  dipole  terms  <(>(x)3G(£,x)/3n  and 
<)>(x)3G(£,x)/9x  in  the  hull  and  waterline  integrals,  respectively,  in  the  classical 
integro-dif f erential  equation  defined  by  Equations  (4.10c),  (4.8),  and  (4.9),  take 
the  forms  [<Kx)-<K£)  ]3G(f,x)/3n  and  [<|>(x)-c{>(£)]3G(£,x)/3x,  respectively,  in  the 

modified  integro-dif f erential  equation  obtained  in  this  study.  These  modified  dipole 

■+ 

terms  are  nonsingular,  i.e.,  remain  finite,  as  the  integration  point  x  approaches 
any  field  point  %  where  the  hull  is  smooth  (i.e.,  has  a  tangent  plane). 

A  recurrence  relation  is  proposed  for  solving  the  integro-dif f erential  Equation 
(5.1)  iteratively.  This  recurrence  relation  is  defined  by  Equation  (5.7),  where  the 
initial  (zeroth)  approximation  is  taken  as  the  nonhomogeneous  term  ip(t)  in  the 
integro-dif f erential  Equation  (5.1).  In  the  particular  case  of  potential  flow  about 
an  ellipsoid  (with  arbitrary  beam-to-length  and  draf t-to-length  ratios)  in  trans- 
latory  motion,  along  any  direction,  in  an  unbounded  fluid,  the  first  approximation, 
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given  by  Equation  (5.6),  actually  is  exact,  as  is  proved  in  Reference  11;  this  first 
approximation  was  also  shown  to  provide  a  good  approximation  to  the  exact  potential 
for  arbitrary  translatory  motions  in  an  unbounded  fluid  of  a  cylinder  in  the  shape 
of  an  ogive  with  arbitrary  thickness  ratio. 

The  plan  of  the  study  is  as  follows.  The  basic  potential-flow  problem  of  the 
three-dimensional  theory  of  flow  about  a  ship  moving  with  constant  speed  in  regular 
waves  is  briefly  formulated  in  Section  2;  a  more  detailed  formulation  of  the  problem 
may  be  found  elsewhere,  for  instance  in  Reference  1.  The  basic  equations  satisfied 
by  the  Green  function  associated  with  the  free-surface  boundary  condition  (2.5)  are 
given  in  Section  3.  Specifically,  the  Green  function,  G,  satisfies  Equations  (3.3a 
and  b)  or  (3.4a  and  b) ,  depending  upon  whether  the  singularity  is  fully  submerged 
(j,<0)  or  exactly  at  the  mean  sea  surface  (?=0),  respectively.  Equations  (3.3a  and  b) 
for  a  fully  submerged  source  are  well  known.  However,  Equations  (3.4a  and  b),  cor¬ 
responding  to  a  flux  across  the  mean  sea  surface,  are  proper  in  the  limiting  case 
when  the  singularity  is  exactly  at  the  mean  sea  surface.  Equations  (3.3a  and  b)  and 

(3.4a  and  b)  generalize  the  corresponding  equations  obtained  previously  for  the  par- 

12 

ticular  cases  of  ship  wave  resistance  and  of  wave  radiation  and  diffraction  at 
zero  forward  speed.13  Equations  (3.3a  and  b)  and  (3.4a  and  b)  are  used  in  Section  4 
for  obtaining  basic  integral  identities  satisfied  by  the  velocity  potential.  The 
three  classical  identities  (4.10a,  b,  and  c)— valid  strictly  outside,  inside,  and  on 
the  ship  surface,  respectively — are  obtained  first.  However,  the  main  new  result  of 
Section  4  is  identity  (4.13).  This  identity  is  valid  outside,  inside,  and  exactly  on 
the  hull  surface,  and  indeed  is  equivalent  to  the  set  of  the  three  usual  identities 
(4.10a,  b,  and  c) .  The  integral  identity  (4.13)  yields  an  integro-dif f erential 
equation  for  determining  the  velocity  potential  on  the  surface  of  a  ship  moving  at 
constant  speed  in  regular  waves.  This  equation  is  examined  in  Section  5.  Finally, 
an  approach  to  the  numerical  evaluation  of  the  iterative  approximations  defined  in 
Section  5  is  presented  in  Section  6. 
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2.  THE  BASIC  POTENTIAL- FLOW  PROBLEM 

The  basic  potential-flow  problem  of  the  linearized  theory  of  ship  motions  in  a 
regular  sea  is  briefly  formulated  in  this  section.  The  sea  is  assumed  to  be  of 
infinite  depth  and  horizontal  extent.  Water  is  regarded  as  homogeneous  and  incom¬ 
pressible,  with  density  p.  Viscosity  effects  are  ignored,  and  irrotational  flow  is 
assumed.  Surface  tension,  wavebreaking,  spray  formation  at  the  ship  bow,  and  non- 
linearities  in  the  sea-surface  boundary  condition  are  neglected.  A  moving  system  of 
coordinates  (X,Y,Z)  in  steady  translation  with  the  mean  forward  velocity  U  of  the 
ship  is  defined.  Specifically,  the  mean  (undisturbed)  sea  surface  and  the  center- 
plane  of  the  ship  in  its  mean  position  are  taken  as  the  planes  Z  =  0  and  Y  =  0, 
respectively;  the  Z  axis  is  directed  vertically  upwards,  and  the  X  axis  is  directed 
toward  the  ship  bow. 

In  the  above-defined  translating  system  of  coordinates,  the  linearized  sea- 
surface  boundary  condition  takes  the  form 

[g3z+(U9x-8T)2]<r  =  (U9X-3T)P'/P  -  gQ'  on  Z  =  0  (2.1) 

where  g  is  the  acceleration  of  gravity,  T  is  the  time,  $  =  $  (X,T)  is  the  velocity 

potential,  P"  =  P'(X,Y,T)  and  Q"  =  Q'(X,Y,T)  correspond  to  distributions  of  pressure 
and  flux,  respectively,  at  the  sea  surface  (we  have  Q'  =  0  for  all  practical  applies' 
tions,  and  ?'  =  0  except  for  surface-effect  ships),  and  the  notation  3.,,  3X>  3^,  is 
meant  for  the  differential  operators  3/3Z,  3/3X,  3/3T,  respectively. 

The  present  study  is  concerned  with  flows  simple-harmonic  in  time,  with  radiant 

f.  aquency  (o  where  (o  is  the  frequency  of  encounter.  However,  such  free-surface  flows 

are  not  completely  (or  uniquely)  determined  unless  an  appropriate  "radiation  condi- 

14 

tion"  is  imposed,  as  is  well  known  and  is  discussed  by  Stoker,  for  instance.  A 
convenient  alternative  approach,  employed  previously  in  Lighthill  and 
Noblesse,12,13  to  the  use  of  such  a  "radiation  condition"  consists  in  defining  a 
time-harmonic  flow  as  the  limit,  as  the  small  positive  auxiliary  parameter  a 
vanishes,  of  a  flow  defined  by  a  velocity  potential  of  the  form 

4>'(X,T)  =  Re<t>  (X)exp  [  (a-ico)T]  (2.2) 
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where  Re  represents  the  real  part  of  the  function  on  the  right  side.  The  eventual 
sea-surface  distribution  of  pressure  P'CX.Y.T)  and  flux  (^(X.Y.T)  similarly  are 
assumed  to  be  of  the  form 

P'(X,Y,T)  =  Re  P(X,Y)exp[(o-iaj)T]  (2.2a) 

Q'(X,Y,T)  =  Re  Q(X,Y) exp [ (a-iw)T]  (2.2b) 

In  this  alternative  approach,  one  is  then  faced  with  an  initial-value  problem,  with 
the  obvious  initial  conditions  =  0  and  3$>V3T  =  0  for  T  =  -“.  Use  of  Equations 

(2.2)  and  (2.2a  and  b)  in  Equation  (2.1)  then  yields  the  sea-surface  boundary 
condition 


[g8z-(aj-iU3x+ia)2]<I>  =  i(a)-iU3x+ia)P/p  -  gQ  on  Z  =  0 


(2.3) 


for  the  "spatial  component"  $(X)  of  the  actual  potential  3>'(X,T). 

Nondimensional  variables  are  defined  in  terms  of  l/o)  as  reference  time,  the 

ship  length  L  as  reference  length,  and  the  acceleration  of  gravity  g  as  reference 

1/2  1/2 

acceleration,  from  which  the  reference  velocity  (gL)  ,  potential  (gL)  L,  and 
pressure  pgL  can  be  formed.  The  nondimensional  variables 

t  =  ft)T,  x  =  X/L,  $  =  <t>/(gL)1/2L,  p  =  P/ pgL,  q  =  Q/(gL)1/2  (2.4) 

are  then  defined.  In  terms  of  these  nondimensional  variables,  the  sea-surface 
boundary  condition  (2.3)  can  be  shown  to  take  the  fonn 

[3z-(f-iF3x+ie)2]<j>  =  i(f-iF3x+i£)p  -  q  on  z  =  0  (2.5) 


where  f  is  the  frequency  parameter,  F  is  the  Froude  number,  and  £  is  the  time- 
growth  parameter  defined  as 


f  =  0)(L/g) 


(2.6a) 


F  =  U/(gL)J 


(2.6b) 


e  =  0 (L/ g) 


(2.6c) 


The  basic  potential-flow  problem  of  the  linearized  theory  of  ship  motions  in  a 
regular  sea  may  now  be  stated.  As  is  well  known,  the  problem  consists  in  solving 
the  Laplace  equation 

V2<J>  =  0  in  d  (2.7) 

subject  to  the  boundary  conditions  specified  below.  The  solution  domain  d  in  Equa¬ 
tion  (2.7)  is  the  domain  exterior  to  the  ship  hull  and  bounded  upwards  by  the  mean 
sea  surface  0.  On  the  mean  sea  surface  a,  the  boundary  condition  (2.5)  must  be 
satisfied: 


[9  -(f-iF9x+ie)  =  i(f-iF9x+ie)p  -  q  on  a 


(2.8) 


where  we  generally  have  p  =  0  =  q.  The  potential  <f>(x)  vanishes  as  |  |x|  |-*  °°  at  least 
as  fast  as  l/||x||;  that  is,  we  have 

4>=0(l/|[x||)as||x||->-°o  (2.9) 

Finally,  on  the  mean  position  of  the  ship  hull  surface  h  the  potential  must  satisfy 
the  usual  Neumann  condition 


9<}>/9n  given  on  h 


(2.10) 


where  94>/9n  =  V<J>*n  is  the  derivative  of  <J>  in  the  direction  of  the  unit  normal  vector 
n  to  h,  taken  to  be  pointing  inside  the  fluid.  The  precise  form  taken  by  the 
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expression  for  34>/3n  on  h  in  the  usual  "radiation"  and  "diffraction"  problems  may  be 
found  in  Newman,'*'  for  instance. 

A  classical  technique  for  solving  a  potential-flow  problem  such  as  that  defined 
by  Equations  (2.7)  through  (2.10),  in  the  general  case  of  an  arbitrary  ship  form, 
consists  in  formulating  an  integral  equation  for  the  potential  based  on  the  use  of 
a  Green  function  satisfying  all  the  equations  of  the  problem  except  the  "hull  boun¬ 
dary  condition,"  which  is  to  be  satisfied  by  means  of  the  integral  equation.  The 
required  Green  function  is  defined  in  the  following  section. 
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3.  THE  GREEN  FUNCTION 

“►  —► 

The  Green  function,  G(£,x),  associated  with  the  sea— surface  boundary  condition 
(2.5)  satisfies  the  equations 

G  =  6(^-x)6(n-y)6(?-z)  in  £  <  0  (3.1a) 

[3^-(f-iF3^+i£)2]G  =  0  on  (  =  0  (3.1b) 

where  6(  )  is  the  usual  Dirac  "delta  function,"  and  represents  the  differential 
operator  (3^, 3^, 3^).  Physically,  the  Green  function  G(£,x)  is  the  "spatial  compo¬ 
nent"  of  the  velocity  potential  ReG(f,x)exp [ (e/f-i) t ]  of  the  flow  created  at  the 
field  point  £(£,n>  ££0)  by  a  moving  source  of  pulsating  strength  Re  exp [ (£/f-i) t ] 
located  at  point  x(x,y,z<0).  In  the  limiting  case,  z  =  0,  the  source  at  point  x 
evidently  is  no  longer  fully  submerged,  so  that  this  physical  interpretation  of  the 
Green  function  becomes  ambiguous.  A  complementary  physical  interpretation  for  this 
limiting  case  is  that  the  pulsating  flow  created  at  point  x(x,y,z=0)  stems  from  a 
flux  across  the  plane  z  =  0  of  the  mean  sea  surface.  In  the  limit  z  =  0,  the  Green 
function  G(£,x)  must  then  satisfy  the  equations 

V2G  =  0  in  £  <  0  (3.2a) 

[3^-(f-iF3^+ie) 2 ]G  =  -6 (£-x)6 (n-y)  on  £  =  0  (3.2b) 

as  may  be  seen  from  the  sea-surface  boundary  condition  (2.5).  Equations  (3.2a  and 

b),  justified  above  on  physical  grounds,  can  be  justified  mathematically  in  the  man- 

12  13 

ner  shown  in  Noblesse  ’  for  the  particular  problems  of  wave  radiation  and  diffrac¬ 
tion  at  zero  mean  forward  speed  (F=0)  and  of  steady  flow  about  a  ship  advancing  at 
constant  speed  in  calm  water  (f=0) . 

The  Green  function  G(^,x)  actually  is  a  function  of  the  four  variables  £-x, 

2 

p-y,  £+z,  and  (£-z)  ,  and  thus  is  invariant  under  the  substitutions  £  *-*■  -  x, 

H  -  y,  £  z.  By  performing  these  changes  of  variables  in  Equations  (3.1a  and  b) 
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and  (3.2a  and  b) ,  it  then  may  be  seen  that  the  Green  function  G(f,x)  also  satisf: 
the  following  equations 


V  G  =  6(x-£)<5(y-n)6(z-0  in  z  <  0 


(3.3a) 


for  C  <  0 


(3  -(f+iF3  +ie)  ]G  =  0  on  z  =  0 
z  x 


(3.3b) 


V  G  =  0  in  z  <  0 


[ 3z~ (f+iFS^+ie)  ]G  =  -<$(x-£)5(y-ri)  on  z  =  0 


for  C,  =  0 


(3.4a) 


(3.4b) 


where  V  is  the  differential  operator  (3^,3  ,3z).  Equations  (3.3a  and  b)  and  (3.4a 
and  b)  will  be  used  in  the  next  section  for  obtaining  integral  identities  satisfied 
by  the  velocity  potential. 

A  well,  known  expression  for  the  Green  function,  in  terms  of  a  double  integral, 
can  be  obtained  by  using  a  double  Fourier  transformation  of  Equations  (3.1a  and  b) 
with  respect  to  the  horizontal  coordinates  £  and  ri>  This  Fourier  representation  of 
the  Green  function  is  given  by 

-+  -*  2  2  2  -1/2  2  2  2  -1/2 
4ttg(C,x)  =  -[(x-O  +(y-n)  +(z-0  ]  L,i  +  l(x-5r+(y-n)  +(z+«  I  ' 


■\  \ 


exp  [ (z+C)  (y2+v2)  (x-Qp+(y-n)v)  ] 

(y2+v2)1/2  -  (f-Fy+ie)2 


(3.5) 


The  "Cartesian  Fourier  integral  representation"  (3.5)  can  also  be  expressed  in 
the  form  of  a  "polar  Fourier  representation"  by  performing  the  change  of  variables 
y  =  X  cos0  and  v  =  X  sinQ,  which  express  the  Cartesian  Fourier  variable  y  and  v  in 
terms  of  the  polar  variables  X  and  0.  These  equivalent  double-integral  representa¬ 
tions  were  first  obtained  by  Haskind^  and  Brard,^  and  later  by  Hanaoka, 
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Stretenski,  Eggers,  Havelock,  and  Wehausen,  and  are  therefore  well  known. 
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More  recently,  one-fold  Integral  representations  (involving  the  exponential 

23 

integral  in  the  integrand)  have  been  obtained  and  used  by  Inglis  and  Price  and 
3 

Guevel  and  Bougis.  These  single- integral  representations  are  modifications  of 
the  double- integral  Fourier  representation  in  terms  of  the  polar  coordinates  (A,0). 
Single- integral  representations  associated  with  the  Cartesian  Fourier  representation 
(3.5),  in  the  manner  shown  in  Reference  12  for  the  particular  problem  of  ship  wave 
resistance  (f=0,  F^O),  have  not  been  obtained  to  the  author's  knowledge.  However, 
such  single- integral  representations  are  considerably  more  complex  than  the  corre¬ 
sponding  integral  representations  for  the  ship  wave  resistance  problem  and  the  series 
representations  obtained  in  Reference  13  for  the  particular  case  of  wave  radiation 
and  diffraction  at  zero  forward  speed  (F=0,  f^O).  For  the  practical  purpose  of 
numerically  evaluating  the  velocity  potential  defined  by  a  surface  (or  line)  distri¬ 
bution  of  singularities  (sources  or  dipoles)  with  known  strength,  it  may  actually  be 
preferable  to  use  a  double- integral  Fourier  representation,  such  as  that  given  by 
Equation  (3.5),  together  with  an  interchange  in  the  order  of  integration  between  the 
Fourier  variables  (p,v)  and  the  space  variables  (x,y,z),  as  is  shown  in  Section  6. 
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4.  FUNDAMENTAL  INTEGRAL  IDENTITIES 

In  this  section,  basic  integral  identities  for  the  velocity  potential  are 
obtained  by  applying  a  classical  Green  identity  to  the  potential  <J>  5  $(x)  and  the 
previously  defined  Green  function  G  5  G(f,x).  The  Green  identity  is 


I 


(4>V2G-GV2<J))dv  = 


i 


(4>3G/3z-G3<})/3z)dxdy 


+ 


|  (G3i}>/3n-<j)3G/3n)da  + 
h 


J  (<J>3G/3n-G3<J>/3n)da 


(4.1) 


where  d'  is  the  finite  domain  bounded  by  the  ship  hull  surface  h,  the  mean  sea  plane 
z  =  0,  and  some  arbitrary,  but  sufficiently  large,  exterior  surface  h^  surrounding 
the  ship  surface  h,  as  is  shown  in  Figure  1;  furthermore,  o'  is  the  portion  of  the 
plane  z  =  0  between  the  intersection  curves  c  and  c^  of  the  plane  z  =  0  with  the 
ship  surface  h  and  the  exterior  surface  h^,  respectively.  On  the  surfaces  h  and  h 
we  have  3<}>/3n  =  V<|>*n  and  3G/3n  =  VG*n  where  n  is  the  unit  outward  normal  vector  to 
h  or  h^,  as  is  shown  in  Figure  1.  Finally,  dv  and  da  represent  the  differential 
elements  of  volume  and  area  at  the  integration  point  x  of  the  domain  d'  and  the  sur¬ 
faces  h  or  h  ,  respectively,  and  dxdy  is  the  differential  element  of  area  of  the 
mean  sea  surface  o'. 

Let  the  integrand  <p 3G/3z  -  G3<J>/3z  of  the  sea-surface  integral  in  Equation  (4.1) 
be  expressed  in  the  form  <t>[9z~(f+iF3x+i£)  2]G  -  G[3z~(f-iF3x+ie)2](}>  + 
2iF(f+ic)3(G<J>)/3x  +  F23(G3<J>/3x-<J)3G/3x)/3x.  Furthermore,  we  may  use  the  relation 
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Figure  1  -  Definition  Sketch 


/ 


[2i(f+ie)9(G<t>)/9x+F9(G9<}>/9x-<|>9G/9x)/9x]dxdy 


[2i(f+ie)G4H-F(G9({>/9x-4>9G/9x)  ]dy 


+ 


[2i(f+ie)G<t>+F(G9<|>/9x-<{>9G/9x)  ]dy 


(4.2) 


where  the  curves  c  and  c^  are  oriented  clockwise  and  counterclockwise,  respectively, 
as  is  shown  in  Figure  1.  The  Green  identity  (4.1)  can  then  be  expressed  in  the  form 
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f  4>V2Gdv  -  f  <}>  [  3  z~  (f +iF3x+ic ) 2  ] Gdxdy 

a 


-  j*  GV  <J)dv  -  j  G [ 3^- (f -iF3^+i£) 2 ] $dxdy 

j  - 


+  (G3<f>/3n-<{>3G/3n)da 


+  F  J  [2i(f+ic)G<}>+F(G3<<>/3x-$3G/3x)]dy  +  Io 


(4.3) 


where  the  term  IM  is  given  by  the  integrals 


!„  =  J  (<pdG/dn-Gd<p/dn)da  +  F  j"  [2i(fidx)G<M-F(G3<t>/3x-c|>3G/3x)  ]dy 


7  ?  7  l  /  7 

We  have  G  =  0(l/r)  and  <J>  =  0(l/r)  as  r  =  (C  +n  +?  )  -*■  so  that  the  term  1^ 

vanishes  as  the  large  surrounding  surface  h^  is  made  ever  larger.  The  term  1^  can 
then  be  ignored  if  the  finite  domain  d'  and  the  finite  region  a'  of  the  mean  sea 
plane  are  replaced  by  the  unbounded  mean  flow  domain  d  and  the  unbounded  mean  sea 
s'.irfa^o  n  outside  the  mean  hull  surface  h  and  its  intersection  curve  c  with  the  plane 
z  =  0,  respectively. 

By  expressing  the  potential  <f>  in  the  integrands  of  the  two  integrals  on  the  left 
side  of  Equation  (4.3)  in  the  form  <}>  =  +  (4|-41*) ,  where  <p  =  <p(x)  as  was  defined 

previously,  and  <t>*  represents  the  potential  at  the  field  point  t,,  i.e.,  =  <K£)> 

we  may  obtain 


|  4>V2Gdv  -  j"  4>[32-(f+iF3x+ic)2]Gdxdy  =  C<}>*  +  C' 

d  a 


I 


where  C  and  C'  are  defined  as 


C  =  j  V2Gdv  -  j  [3z-(f+iF9x+ie)2]Gdxdy 
d  a 

C'  =  J  V2G dv  -J  <0-(f>*)  [32-(f+iF3x+ic)2]Gdxdy  (4.5) 


It  may  be  seen  from  Equations  (3.3)  and  (3.4)  that  we  have  C'  =  0  if  cJj  -  <J)  =  <J>(x) 

-  <t>U)  +  0  as  x  +  that  is  if  the  potential  is  continuous  everywhere  in  the  solu¬ 
tion  domain  d  and  on  its  boundary  a  +  h  +  c,  as  is  assumed  here.  Use  of  Equation 
(4.4),  with  C'  =  0,  in  Equation  (4.3)  then  yields 


C<)>^  =  J  GV2c(>dv  -  j  G[3z-(f-iF9x+ie)2]0dxdy 
d  a 

+  J  (G3<|>/3n-<J>3G/3n)da 
h 


+  F 


J 


[2i(f+ie)G<M-F(G34>/3x-4>3G/3x)  ]dy 


(4.6) 


Let  t(tx>t^,0)  represent  the  unit  vector  tangent  to  the  curve  c  oriented  in  the 

clockwise  direction,  as  is  shown  in  Figure  1.  On  the  mean  waterline  c,  we  have 

dy  =  t  dH,  where  d £  is  the  differential  element  of  arc  length  of  c.  Furthermore,  we 

have  3<$>/3x  =  V<}>*i,  where  i  (1,0,0)  is  the  unit  positive  vector  along  the  x  axis.  This 

then  yields  3(J>/3x  =  [n3<j>/3n+t3<t>/3£+(n><t)3<)>/3d]  •  i  =  n  3 4>/ 3n  +  t  3<j>/3&  -  n  t  3<(>/3d, 

x  x  z  y 

where  (n  ,n  »nz)  are  the  components  of  the  unit  outward  normal  vector  n  to  the  hull 
surface  h,  34>/3 Jl  is  the  derivative  of  4>  in  the  direction  of  the  tangent  vector  t  to 

-V  •> 

c,  and  3<J)/3d  is  the  derivative  of  <j>  in  the  direction  of  the  unit  vector  n  *  t,  which 
is  tangent  to  h  and  pointing  downwards. 
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Equation  (4.6)  can  then  be  expressed  in  the  form 

C<p (t)  =  I Hi)  -  Ut-A)  (4.7) 

where  C  is  given  by  Equation  (4.5),  ijj(|)  is  the  potential  defined  as 

'HO  =  j  GV  <})dv  -  j  G  [  3^- (f-iFS^+ie) 2  ]  <})dxdy 
d  a 

+  J"  G3(J)/3nda  +  F2  J  Gnxty3<j>/3nd£  (4.8; 

h  c 

and  L(£;$)  is  the  linear  transform  of  tp  defined  as 

L(f;<f>)  =  J  4>3G/3nda  -  2i(f+ie)F  J  C<ptydZ 

h  c 

+  F2  [$3G/3x-G(t  34>/3S,-n  t  3$/3d)]t  d£  (4.9) 

J  x  z  y  y 

c 

Use  of  Equations  (3.3a  and  b)  and  (3.4a  and  b)  in  expression  (4.5)  for  C  then  shows 
_  ■+ 

that  we  have  C  -  1  if  the  field  point  £  is  strictly  outside  the  hull  surface  h,  in  d 

_ 

or  on  o,  whereas  we  have  C  =  0  if  £  is  strictly  inside  the  ship  surface  h.  It  can 
also  be  seen  from  Equations  (3.3)  and  (3.4)  that  we  have  C  =  1/2  if  the  point  £  is 
exactly  on  the  hull  surface  h  or  on  its  intersection  c  with  the  plane  z  =  0,  at 
least  for  points  t  where  the  hull  h  +  c  is  smooth;  more  generally,  the  value  of 
4ttC  (or  2ttC)  at  a  point  t  of  h  (or  c)  is  equal  to  the  angle  at  which  d  (or  a)  is 
viewed  from  the  point  %.  We  thus  have 
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'  ' 

' 

.  ."t.  .■>  -> 

ind+o-h-c 

(4.10a) 

0 

"  =  >K£)  -  L(5;4>)  for  £  < 

in  dj  +  0^  -  h  -  c 

►  (4.10b) 

(p(t)f  2 

exactly  on  h  +  c 

,  d 

(4.10c) 

where  d^  and  cj^  represent  the  domain  and  the  portion  of  the  plane  z  -  0,  respec¬ 
tively,  strictly  inside  the  ship  surface  h,  as  is  shown  in  Figure  1. 

The  value  of  the  constant  C  on  the  left  side  of  Equation  (A. 7)  is  discontinuous 
across  the  ship  hull  surface  h;  C  being  equal  to  1  outside  h  and  to  0  inside,  as  is 
explicitly  indicated  in  Equations  (4.10a,  b,  and  c).  This  discontinuity  in  the  value 
of  C  evidently  is  accompanied  by  a  corresponding  discontinuity  on  the  right  side  of 
Equation  (4.7).  Specifically,  the  latter  discontinuity  stems  from  the  dipole- 
distribution  integrals  <}>3G/3nda  and  4)3G/9xt^di!.  in  the  potential  L(£;<j>) 

defined  by  Equation  (4.9).  An  identity  valid  for  any  point  K — outside,  inside,  or 
exactly  on  the  ship  surface  h — can  be  obtained  by  eliminating  the  discontinuity  in 
the  value  of  C  in  Equation  (4.7).  This  can  be  done  by  adding  the  term  on  both 

sides  of  Equation  (4.7),  with  given  by 


C 


i 


V2Gdv  - 


[  8  - ( f ■ +iF3x+ie ) 2 ] Gdxdy 


(4.11) 


Use  of  the  divergence  theorem 


I 


i 


V2Gdv 


I 


9G/3zdxdy  + 


J 


3G/3nda 
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yields 
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and  is  the  linear  transform  of  <J>  defined  as 


/ 


(<M>*)3G/3nda  -  2i(f+ie)F 


J  G ((M*)  t  dt 

c 


+  F2  I  [(<M>*)3G/3x-G(t  S(t>/an-n  t  3<J>/3d)]t  d«.  (4.15) 

J  *  x  zy  y 

c 

in  which  we  have  <J)  -  (x )  and  4>*  =  <J>(£)  as  was  defined  previously.  In  obtaining 
Equation  (4.13),  the  relation  C  +  =  1  was  used.  This  relation  can  be  obtained  by 

using  Equations  (3.3a  and  b)  and  (3.4a  and  b)  in  Equations  (4.5)  and  (4.11),  which 
yield 


C  +  C  - 


-  f  V2Gdv  -  I* 
z<0  z=0 


[3z-(f+iF3x+ie)  ]Gdxdy 


Identity  (4.13)  is  valid  for  any  point  whether  outside,  inside,  or  exactly  on  the 

ship  surface  h.  This  identity  thus  is  essentially  equivalent  to  the  set  of  the  three 

classical  identities  (4.10a,  b,  and  c),  which  are  exclusively  valid  for  £  outside, 

inside,  and  on  the  hull  surface  h,  respectively. 

Identities  (4.10a,  b,  and  c)  and  (4.13)  correspond  to  the  case  of  an  open  hull 

surface  piercing  the  sea  surface.  For  a  closed,  fully  submerged  surface  h,  the 

waterplane  integral  w(£)  defined  by  Equation  (4.14)  and  the  integrals  around  the 

mean  waterline  c  in  Equation  (4.8)  for  the  potential  ^(f)  and  in  Equations  (4.9) 

and  (4,15)  for  the  potentials  L(£;4>)  and  L'(f;<|>)  are  evidently  not  present.  Two 

other  important  particular  cases  of  identities  (4.10a,  b,  and  c)  and  (4.13)  are 

obtained  in  the  limiting  cases  when  the  Froude  number  F  vanishes,  corresponding  to 

9 

wave  radiation  and  diffraction  by  a  body  with  zero  mean  forward  speed,  and  when 
the  frequency  parameter  f  vanishes,  corresponding  to  steady  flow  about  a  ship 
advancing  in  calm  water. ^ 
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5.  INTEGRAL  EQUATION  AND  RELATED  ITERATIVE  APPROXIMATIONS 
Identities  (4.10a,  b,  and  c)  and  (4.13)  hold  for  any  function  <j>  continuous  in 
the  domain  d  and  on  its  boundary  a  +  h  +  c.  If  the  function  <|>  is  taken  as  the  veloc¬ 
ity  potential  of  flow  about  a  ship  advancing  at  constant  mean  speed  in  a  regular  sea, 
then  the  normal  derivative  3<t>/3n  of  4>  is  given  on  the  hull  surface  h,  and  <)>  satisfies 
the  Laplace  Equation  (2.7)  in  the  mean  flow  domain  d  and  the  sea-surface  boundary 
condition  (2.8),  with  p  =  0  =  q,  at  the  mean  sea  surface  o.  Identities  (4.10c)  and 
(4.13)  then  yield  integro-dif ferential  equations  for  determining  the  potential  <t>  on 
the  ship  surface.  Specifically,  Equation  (4.13)  becomes 

[1-w(|)]<k|)  =  <P(&  -  L '(£;♦)  (5.1) 

where  the  waterplane  integral  w(£)  is  given  by 


+  F2  f[(<M>*)8G/9x-G(t  3*/3t-n  t  3<*/3d)]t  d* 
J  x  z  y  y 

c 


(5.4) 
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The  potential  4<(C)  is  known  since  3 4>/ 3n  is  given  on  the  ship  surface  h.  The  poten¬ 
tial  1/ (!;$),  on  the  other  hand,  is  evidently  not  known. 

An  approximate  solution  of  the  integro-dif f erential  Equation  (5.1)  may  be 
obtained  by  seeking  a  solution  of  Equation  (5.1)  of  the  form  <}>(f)  =  k(f)iKf),  where 
the  function  k(f)  =  <K^)/’K^)  is  assumed  to  be  slowly  varying.  Specifically,  by 
adding  the  term  k(Ql/  (£;i|j)  to  both  sides  of  Equation  (5.1)  and  multiplying  the 
resulting  equation  by  <p(0,  we  may  obtain 


<Kf)Ul-w(f)]<Kl)+l/(f;ijj)}  =  <^2(l)  +  <J>(!)l/(f;i|/)  -  (5.5) 

If  the  potential  <}>  were  actually  proportional  to  the  potential  \Jj,  the  term 
(O L  (£;ijj)  -  <K£)L  (4;<t>)  would  vanish,  and  the  modified  integro-dif f erential  Equa¬ 
tion  (5.5)  would  yield  the  solution 

<P(t)  “  <p2(t)/([l-w(t)l>p(t)  +  L'(t;ip>}  (5.6) 

More  generally,  the  above  expression  for  the  potential  can  be  regarded  as  the 
first  approximation  in  the  sequence  of  iterative  approximations  associated 

with  the  recurrence  relation 

<t>(n+1)(t)  =  <Kl)<t><n)(l)/Ul-w(l)]<t>(n)(l)  +  I/(f;<!>(n))}  for  n  >  0  (5.7) 

and  the  initial  (zeroth)  approximation  4>^(£)  =  ty(C)  •  An  approach  to  the  numerical 
evaluation  of  these  Iterative  approximations  is  presented  in  the  following  section. 
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6.  NUMERICAL  EVALUATION  OF  ITERATIVE  APPROXIMATIONS 
Equation  (3.5)  for  the  Green  function  may  be  written  in  the  form 


00 

4ttC(|,x)  =  -(1/r-l/O  -  (1/tt)  J  dv  J 


dyE(p,v;£)E(|j,v;x)/D(p,v)  (6.1) 


where  r,  r  ,  E(p,v;x),  E(p,v;0,  and  D(p,v)  are  defined  as 


r  =  [(x-02+(y-n)2+(z-C)2]1/2 


(6.1a) 


r "  =  [(x-?)2+(y-n)2+(z+02]1/2 


(6.1b) 


->  2  2  1/2 
E(p,v;x)  =  exp[z(p  +v  )  +i(xp+y\>) 


(6.1c) 


E(p,V;£)  =  exp(^(p2+v2)1/,2-i(Cu+UV) ; 


(6. Id) 


D(p,v)  =  (p2+v2)1/2-(f-Fp+it)2 


■.  b .  1  e ) 


The  potential  i(i(£)  defined  by  Equation  (5.3)  can  then  be  expressed  in  the 


<Kf)  =  +  /(I) 


S  R 

where  the  potentials  tp  and  ip  correspond  to  the  singular  terms  1/r  -  1/r '  and  the 
regular  term  defined  by  the  double  integral,  respectively,  in  Equation  (6.1)  for 

S 

the  Green  function.  Specifically,  the  potential  ^  is  given  by  the  hull-surface 
integral 


22 


^S(f)  =  -(1/4tt>  J  (l/r-l/r')9<t>/9nda 
h 


(6.2a) 


and  the  potential  iJ1  may  be  expressed  in  the  form 


A!)  =  -(1/4ti2) 


OO 

/  du  J 


dME(p,v;C)A(u,v)/D(p,v)  (6.2b) 


where  A(p,v)  is  defined  as 


A(y,v)  =  J  E(p,v;x)3<|>/3nda  +  F  E(p,v;x)nxt  90/9nd£ 
h  c 


(6.2c) 


The  waterplane  integral  w(£)  defined  by  Equation  (5.2)  takes  the  form 


w(f)  = 


(f+i£)' 

i  2 


CO  OO 

j  dv  j  dp  J  exp[i(xu+yv)]dxdy 


(6.3) 


OO  —  <30 


The  potential  L'(£;ip)  defined  by  Equation  (5.4)  can  then  be  expressed  in  the 


l/(f;ij>)  =  Lgd;*)  +  L'(|;i(.) 


(6.4) 


where  the  potentials  L'  and  L'  are  defined  as 

o  K 
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L'(!;40  =  -(1/ 4tt)  [(^S+/)-(^*)]3(l/r-l/r')/8nda 


(6.4a) 


L"(C;40  =  -(l/4ir  )  dv  duE(M,v;OA"(M,v)/D(p,v) 


(6.4b) 


with  A^(y,v))  given  by 


■(U,V)  =  J  [dS+/)-(^*)]3E(h,v;x)/3nda 


2i(f+ie)F  I  (/-^)E(y,V;x)tydA 


+  F2  {dR-^)3E/9x-[t  dipR/dl-n  t  3dS+/)/3d]E}t  d£  (6.4c) 
J  x  z  y  y 


The  first  iterative  approximation  ;(C)  can  then  be  determined  by  using  Equa¬ 
tions  (6.2),  (6.3),  and  (6.4)  in  Equation  (5.6).  The  potential  L'’(C,<f>^n^) ,  n  _>  1, 
defined  by  Equation  (5.7),  for  the  second  and  subsequent  iterative  approximations 
<f>  (n+1)  (n_>l) ,  can  be  expressed  in  a  form  almost  identical  to  that  given  above  by 
Equations  (6.4)  and  (6.4a,  b,  and  c)  for  the  potential  L'(C;4'). 

The  basic  computational  task  common  to  Equations  (6.2b  and  c) ,  (6.3),  and 
(6.4b  and  c)  consists  in  evaluating  a  double  Fourier  integral,  1(C)  say,  of 
the  form 


=«*»  ■  I 


dv  dyexpU(u2+v2)1/2-i(Cy+nv)]N(y,v)/D(M,v)  (6.5) 


-  CO  -  CO 
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where  D(y,v)  is  given  by  Equation  (6.1e),  and  N(y,v)  is  defined  by  a  surface  or  a 
line  integral  of  the  type 


N(p,v) 


■I 


2  2  1/2  -► 
exp[z(u  +v  )  +i(xy+yv)  ]A(x)da,dxdy ,dS. 


(6.5a) 


h,ai,c 


If  a  small  positive  number  is  used  for  the  parameter  e,  the  function  D(y,v)  has  no 
real  root,  and  the  integral  (6.5)  can  be  evaluated  without  difficulty  in  principle. 
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